Bonnor-type Black Dihole Solution 
in Brans-Dicke-Maxwell Theory 

Hongsu Kim* and Hyung Mok Lee^ 
Astronomy Program, SEES, Seoul National University, Seoul, 151-742, KOREA 

Abstract 

It was originally thought that Bonnor's solution in Einstein-Maxwell theory describes a singular 
point-like magnetic dipole. Lately, however, it has been demonstrated that indeed it may describe 
a black dihole, i.e., a pair of static, oppositely-charged extremal black holes with regular horizons. 
Motivated particularly by this new interpretation, in the present work, the construction and exten- 
sive analysis of a solution in the context of the Brans-Dicke-Maxwell theory representing a black 
dihole are attempted. It has been known for some time that the solution-generating algorithm 
of Singh and Rai produces stationary, axisymmetric, charged solutions in Brans-Dicke-Maxwell 
theory from the known such solutions in Einstein-Maxwell theory. Thus this algorithm of Singh 
and Rai's is employed in order to construct a Bonnor-type magnetic black dihole solution in Brans- 
Dicke-Maxwell theory from the known Bonnor solution in Einstein-Maxwell theory. The peculiar 
features of the new solution including internal infinity nature of the symmetry axis and its stability 
issue have been discussed in full detail. 
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I. Introduction 



Not surprisingly, exact solutions of general relativity describing multiple black hole con- 
figurations are few since they would be highly non-trivial to find or construct. Although one 
would naturally expect that such configurations might in general involve a very complicated 
structure, there actually exist some simple solutions having remarkable properties. These 
include the Majumdar-Papapetrou solutions Q| representing an arbitrary number of static, 
extremal charged black holes, all with charges of the same sign and the multi-Schwarzschild 
solution of Israel and Khan Besides, the C and the Ernst metrics ^ or the cosmo- 
logical multi-black hole solutions of Kastor and Trachen P| , all describing the black holes in 
relative motion, may also fall into this category of multi-black hole solutions. 
In the present work, we are particularly interested in Bonnor's solution. Several years ago, 
using a technique to generate a static, axisymmetric solution of the Einstein-Maxwell the- 
ory from the stationary, axisymmetric Kerr solution of the vacuum Einstein theory, Bonnor 
^ constructed a solution describing a magnetic dipole. Bonnor's solution was originally 
thought to describe a singular point-like dipole. Recently, the generalization of Bonnor's 
solution in Einstein-Maxwell theory to its counterparts in Einstein-Maxwell-dilaton theories 

with the correct analysis of the conical singularity structure of the solutions was made by 

n 

Davidson and Gedalin '7']. Interestingly, the new identifications of these solutions with a 
pair of oppositely-charged extremal black holes (rather than with a conventional singular 
point-like dipole) have appeared only afterwards. To our knowledge. Gross and Perry 
were the first to find the Kaluza-Klein analog of the Bonnor's solution and identify it with 
a "self-gravitating" pole-antipole configuration. Later on, it has been further refined by the 
work of Dowker et al. ^ where a thorough analysis of the structure of the solutions and the 
introduction of the background magnetic field were carried out. Even in these works, the 
interpretation was made on the basis of the particular structure of the Kaluza-Klein theory, 
which cannot be applied to solutions with other values of the dilaton. It is also interesting 
to note that the Kaluza-Klein dipole solution of Gross and Perry has lately been uplifted 
to its counterpart in = 11 M-theory by Sen llOll which, via the M/IIA string duality, can 



be interpreted as the D6 — D6 pair solution [11|, [12] in = 10 IIA supergravity theory. 
That the magnetic dipole solution of Einstein-Maxwell-dilaton theories, including of course 
the Bonnor's solution, actually describes a dihole, i.e., a pair of static, oppositely-charged 
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extremal black holes with regular horizons was first demonstrated by Emparan [13]. Moti- 
vated particularly by this work of Emparan, in the present work, we would like to present the 
construction and extensive analysis of a solution in the context of the Brans-Dicke-Maxwell 
theory representing a pair of static, oppositely-charged extremal black holes. As such, it 
might be relevant to mention the current status of Brans-Dicke (BD) theory as a viable the- 
ory of gravity and describe briefly the strategy to be employed to construct a new solution 
of BD-Maxwell theory from a known solution of Einstein-Maxwell theory. 
Perhaps the Brans-Dicke theory is the most studied and hence the best-known of all 
the alternative theories of classical gravity to Einstein's general relativity, . This theory can 
be thought of as a minimal extension of general relativity designed to properly accomodate 
both Mach's principle and Dirac's large number hypothesis 2. Namely, the theory 
employs the viewpoint in which the Newton's constant G is allowed to vary with space 
and time and can be written in terms of a scalar ("BD scalar") field as G = 1/$. Like 
in Einstein's general relativity, to find the exact solutions of the highly non-linear BD field 
equations is a formidable task. For this reason, algorithms generating exact, new solutions 
from the known solutions of simpler situations either of the BD theory or of the conven- 
tional Einstein gravity have been actively looked for and actually quite a few have been 
found. To the best of our knowledge, methods thus far discovered along this line includes 
those of 
and Rai 

to generate stationary, axisymmetric solutions in vacuum BD theory from the known Kerr 
solution in vacuum Einstein theory and later on Singh and Rai generalized this method 
to the one that generates stationary, axisymmetric, charged solutions in BD-Maxwell theory 
from the known Kerr-Newman (KN) solution 19| in Einstein-Maxwell theory. Thus in the 



anis et al., Buchdahl, Mcintosh, Tupper Tiwari and Nayak |r^, and Singh 
In particular, Tiwari and Nayak proposed an algorithm that allows us 



present work, we shall employ the solution-generating algorithm of Singh and Rai in order 
to construct a Bonnor-type magnetic dipole solution in BD-Maxwell theory from the known 
Bonnor solution in Einstein-Maxwell theory and discuss its unfamiliar nature in full detail. 



II. Solution-generating algorithm of Singh and Rai 



We begin by briefly reviewing the a' 



gorithm proposed first by Tiwari and Nayak 
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and 



generalized later by Singh and Rai |l8|. Consider the BD-Maxwell theory described by the 



3 



action 

where $ is the BD scalar field and uo is the generic parameter of the theory. Extremizing this 
action then with respect to the metric g^p, the BD scalar field and the Maxwell gauge 
field (with the field strength F^,^ = V^iA,^ — VuA^j) yields the classical field equations 
given respectively by 

r — R _ In p — ^^T^Af I o„rpBD 

where (2) 



rpBD _ J_ 

and 



V„V"$ = ,^ ^ rf ^ = 0, V^F^'' = 0, V^F'"' = 



(2u; + 3) ■ 

with the last equation being the Bianchi identity and F^,^ = |e^(^Fa/3. And the Einstein- 
Maxwell theory is the — oo limit of this BD-Maxwell theory. As mentioned in the 



introduction, the BD theory can be thought of as a minima! 



extension of general relativity 



designed to properly accomodate both Mach's principle M arrd D,ra.-s large number hy- 

pothesis Namely, the theory employs the viewpoint in which Newton's constant G is 
allowed to vary with space and time and can be written in terms of a scalar ( "BD scalar" ) 
field as G = 1/$. Note also that in the action and hence in the classical field equations, 
there are no direct interactions between the BD scalar field $ and the ordinary matter, i.e., 
the Maxwell gauge field A^. Indeed this is the essential feature of the BD scalar field $ that 
distinguishes it from "dilaton" fields in other scalar-tensor theories such as Kaluza-Klein 
theories or low-energy effective string theories where the dilaton-matter couplings generi- 
cally occur as a result of dimensional reduction. (Here we would like to stress that we shall 
work in the context of original BD theory format not some conformal transformation of it.) 

n 

As a matter of fact, it is the original spirit of BD theory of gravity in which the BD 
scalar field $ is prescribed to remain strictly massless by forbidding its direct interaction 
with matter fields. Now the algorithm of Tiwari and Nayak, and Singh and Rai goes as 
follows. Let the metric for a stationary, axisymmetric, charged solution to Einstein-Maxwell 



field equations take tlie form 

ds' = -e^u^{dt + WEd(t>f + e2('=-"^-)[(rfxi)2 + [dx^f] + hle-^""^ dcf)'' (3) 

wliile tlie metric for a stationary, axisymmetric, cfiarged solution to BD-Maxwell field equa- 
tions be 

ds^ = -e^^^^idt + WBDd(t)f + e^^''^^-'^^^\{dx^f + [dx^f] + /^l^e-^^^^rf^^ (4) 

where U , W, k and h are functions of x^ and only. The significance of the choice of the 
metric in this form has been thoroughly discussed by Matzner and Misner [2c| and Misra and 
Pandey j21||]. Tiwari and Nayak, and Singh and Rai first wrote down the Einstein-Maxwell 
and BD-Maxwell field equations for the choice of metrics in eq.(3) and (4) respectively. 
Comparing the two sets of field equations closely, they realized that stationary, axisymmetric 
solutions of the BD-Maxwell field equations are obtainable from those of Einstein-Maxwell 
field equations provided certain relations between metric functions hold. 
That is, if {We, kE, Ue, tiE, A^) form a stationary, axisymmetric solution to the Einstein- 
Maxwell field equations for the metric in eq.(3), then a corresponding stationary, axisym- 
metric solution to the BD-Maxwell field equations for the metric in eq.(4) is given by 
{Wbd, ksD, Ubd, hsD, A^^) where 

Wbd = We, kBD = kE, Ubd = Ue - ^\og(!>, (5) 

hBD = $ = [hEf^'''~'^\ Af" = A^ 

In the following subsections, we shall apply this solution-generating algorithm to the 
construction of the Brans-Dicke gravity Minkowski spacetime (i.e., the empty spacetime 
solution of the Brans-Dicke theory) and the Bonnor-type magnetic dipole solution of 
Brans-Dicke-Maxwell theory. 

(1) Empty spacetime solution in Brans-Dicke theory 

We begin with the construction and interpretation of the empty spacetime solution 
of Brans-Dicke theory which will serve as a reference background geometry for all the 
other spacetime solutions of Brans-Dicke-Maxwell theory. Then it will allow us eventually 
to establish the general physical interpretation of all the new solutions that result upon 
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performing the solution-generating algorithm presented above to the known solutions in 
general relativity. Consider the Minkowski spacetime metric which clearly is a solution to 
the vacuum Einstein equation 

ds'' = -de + dr" + rHe'' + sin^ ^^0^ (g) 

For later comparison with its counterpart in Brans-Dicke theory, we also point out that 
the Newton's constant G is indeed a constant in this Einstein theory. In order to set this 
Minkowski metric into the form well prepared for the solution-generating method given 
above, we first consider the change of radial coordinate r — and hence dr — e^dR which 
leads us to rewrite the eq. (6) into 

ds^ = -dt'' + e^^'idR'' + de^) + e^^ sin^ ed(f)''. (7) 

Then identifying it with the standard form given in eq.(3), we can read off the metric 
components as 

e2^^ = 1, We = 0, (8) 
g2fe, ^ e^^, /il^e^^sin^e. 

Now using the solution-generating rule in eq.(5) in the algorithm by Tiwari and Nayak, and 
Singh and Rai, we can now construct the empty spacetime solution of BD theory as 

ds' = -{r'sin'0y'^^'^^'^dt' + {r'sm'ef^'''^'\dr' + r'de'] 

+ (r'sin'ey'^^'^^'\'sm'ed<p', (9) 

where we restored the "bare" Newton's constant G in the expression for the BD scalar field 
solution. Among other things, for a; — > oo, it reduces to the Minkowski spacetime metric 
with ^ — 1/G — constant in eq.(6) as it should since in this limit the Brans-Dicke theory 
goes over to the Einstein's general relativity. 

Clearly, the solution-generating algorithm of Singh and Rai breaks the spherical symmetry 
along the way from the Minkowski spacetime to the Brans-Dicke gravity empty spacetime 
as one can see in eq.(9) above. Close inspection reveals that this happens due to the emer- 
gence of non-trivial BD scalar field which essentially represents the inverse of the spacetime- 
dependent effective Newton's constant as mentioned earlier, i.e., $(r, ^) = G~jf{r,9). In 
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other words, since the inverse of the BD scalar field representing the strength of the grav- 
itational interaction changes from point to point, both the spherical symmetry and the 
asymptotic fiatness of the usual Minkowski spacetime are broken accordingly. 
We now attempt to provide the physical interpretation of the nature of this empty spacetime 
solution of BD theory in a careful manner. This, however, cannot be achieved unless some 
particular finite value of the BD parameter is specified. Therefore we first need to appreciate 
the physical meaning of the dimensionless BD parameter uj. In the BD gravity action given 
in eq.(l), the term ~ a;(VQ,$V"$/$), like other terms in the action, should be finite. Thus 
large lo indicates the regime where the BD scalar field $ remains nearly constant as the 
inverse of the present value of the Newton's constant. Namely, this regime amounts to the 
Einstein's general relativity limit. On the other hand, the small lo indicates the regime in 
which the BD scalar field varies sizably with space and time and hence the theory deviates 
largely from the general relativity. We have already noted that in the limit a; — > oo, the 
BD gravity empty spacetime in eq.(9) does reduce to the Minkowski spacetime metric with 
^ — IjG — constant as it should. Thus next we consider the nature of this spacetime 
solution in the other limit, namely for small to values. We shall particularly pick the value 
of the 10 parameter to be a; = —2 which lies in the range of later interest, —5/2 < o; < —3/2 
(which has particular significance that will be explained in the next subsection where we 
shall discuss the Bonnor-type magnetic dipole solution in Brans-Dicke-Maxwell theory). For 
this value of the BD a; parameter, the BD gravity empty spacetime solution becomes 

/ 2 • 2 /)\2 / 2 -'- / 2 • 2 a\'i 

9tt = -{r sm 9) , Qrr = gee r = . ^ g^^ = (r sm 9) , 

(r^sm 9y 

$(r,^) = ^(r^sin^^)-^ (10) 

Among others, it is interesting to note that from the behavior of the metric function gu, 
one can realize that near the origin r — and the symmetric axis 9 — 0,7r, the infinite 
time delay occurs and hence all the dynamics (if exist) would be frozen there. This is 
reminiscent of what happens at the horizon of a black hole. Next, from the behavior of the 
metric functions grr and ggo, one can notice that both the origin and the symmetry axis 
are infinite proper distance away ! Namely, they turn out to be internal infinities. Again 
this is very reminiscent of the feature of a black hole horizon which is indeed an internal 
infinity to a distant observer. We now attempt to interpret this peculiar geometric structure 
of the BD gravity empty spacetime solution in terms of the spacetime- dependent behavior 
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of the BD scalar field which, as mentioned earlier, serves as the inverse of an effective 
Newton's constant, $(r, ^) = G~jf{r,9). First, as we just pointed out, both the origin and 
the symmetry axis are infinite proper distance away 

/■O rO 1 

/ \/g^dr = / . oo> (11) 

Jr Jr sm 6 

f Vg^ede = f —^je ^ oo (12) 

Je Je rsm 9 

while the asymptotic region r = oo is finite proper distance away 

/•oo rco ^ 

I V97rdr^ I -^-^^dr< 00. (13) 

It seems that this can be attributed to the spacetime-dependent behavior of the effective 
Newton's constant 

^-\r,e) = Geff{r,e) = G{r^sm^e)^. (14) 

Namely, since the Newton's constant is a measure of the strength of the gravitational inter- 
action, this spacetime-dependent behavior of Gg// implies that the gravitational interaction 
vanishes at the origin and along the symmetry axis while it diverges asymptotically as 
r — > oo. To be a little more concrete, the regions such as the origin and the symmetry axis 
are essentially inaccessible as the gravitational interaction is effectively absent and hence no 
test object would get accelerated (or at least be in motion) there. Meanwhile the coordinate 
infinity r = oo can be in a finite proper distance as the gravitational interaction grows as 
~ and eventually diverges at r = oo and thus any test particle there would get accel- 
erated infinitely. Namely, the physical (proper) distance between the two spacetime points 
gets larger as the effective gravitational interaction there gets weaker whereas it gets smaller 
as the effective interaction there gets stronger. In order to support this type of interpretation 
in terms of the spacetime-dependent behavior of the effective Newton's constant, one can 
take some other values of the BD co parameter and see if essentially the same interpretation 
still holds true. For instance if one takes, say, cj = 1/2 or 5/2 which lie outside of the above- 
mentioned range, —5/2 < a; < —3/2, just the other way around is the case. Namely, the 
origin and the symmetry axis are finite proper distance away (indeed instantly accessible) as 
the gravitational interaction strength diverges there while the asymptotic region r = oo is 
infinite proper distance away as the gravitational interaction strength vanishes there. And of 
course in this analysis, the "backreaction" of the test particle to the background spacetime 



that would modify the background spacetime structure itself is neglected. The peculiar fea- 
tures of the BD gravity empty spacetime geometry, particularly the fact that the symmetry 
axis actually is an internal infinity, in turn, leads to the occurrence of rather an embarrass- 
ing conical singularity and we now turn to this issue. Notice that the conical singularity, or 
more specifically, the deficit angle can be defined in the following manner. Take the ratio 
between the proper circumference and the proper radius of a small limiting circle around 
the symmetry axis 6 = 0,71. If this ratio comes out to be exactly 27r, then there is no deficit 
angle. If instead it turns out to be less than 27r, then there exists non-vanishing deficit angle. 
Being guided by this criterion, we now compute the possible deficit angle in the BD gravity 
empty spacetime geometry. 



^(o.tt) = 27r - I — — — — |e=o,7r = 27r (15) 



where A0 is the period of the azimuthal angle coordinate that will be chosen as 27r and 
we used gee = r^/(r^ sin^ 6^)^ and g^^ = (r^sin^6')'^ and hence {gee)~^^'^d^g^^/ dO\e=Q,TT = 0. 
This demonstrates that this seemingly pathological "27r-angle deficit" is indeed a direct 
consequence of our earlier observation that the proper radius of a small limiting circle 
around the symmetry axis diverges, ^/geedO oo as the symmetry axis can be thought 
of as an internal infinity. Put differently, "27r-angle deficit" is not a generic pathology but 
just another manifestation of the emergence of the internal infinity, i.e., the infinitely far 
symmetry axis. 

(2) Bonnor-type magnetic dipole solution in Brans-Dicke-Maxwell theory 

Next, we apply the same method to obtain the Bonnor-type magnetic dipole solu- 
tion in BD-Maxwell theory from the known Bonnor solution [o^ in Einstein-Maxwell theory. 
And to do so, again one needs some preparation which involves casting the Bonnor solution 



given in Boyer-Lindquist coordinates {t,r,6,(f)) in the metric form in eq.(3) 



oj performing 



a coordinate transformation (of r alone) suggested by Misra and Pandey 21[. Namely, we 
start with Bonnor's magnetic di pol e (dihole) solution of Einstein-Maxwell theory written in 



aiDoie 

m 



Boyer-Lindquist coordinates 

2Mr\2 



ds' = 1 



+ (A + (M2 + a2)sin2^^)3^X + ^ 
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+ (l-'^y\sin^edct)\ (16) 
/\ + sm u 

where E = — a? cos^ 9 and A = — 2Mr — with M denoting the ADM mass and 
a representing (roughly) the proper separation between the poles [l3|. Consider now the 
transformation of the radial coordinate 

R (M^ + a^) p 

r = e^ + M+ -e-^ (17) 

4 

which gives dr^ I IS. = dB?. Indeed, this transformation has been deduced from the associated 

n 

one introduced by Misra and Pandey |21| for the case of single rotating (Kerr) black hole. 
Then Bonnor's magnetic dipole solution can now be cast in the form in eq.(3), i.e., 

E y (A + (M2 + a?) sin^ OfK E 7 ^ ' 

+ (^i_^^"'Asin2^#2 (^g) 

with now E = — cos^ 9 and A = — 2ML — where we set, as a short-hand notation, 
L = + M + (^^ +" ^ Q-R^ Then identifying it with the standard form given in eq.(3), we 
can read off the metric components as 

e-. . (l - H|^)^ (^±i^)\ W..,, (19) 

Now using the solution-generating rule in eq.(5) in the algorithm by Tiwari and Nayak, 
and Singh and Rai, we can now construct the Bonnor-type magnetic dipole solution in 
BD-Maxwell theory as 

2MLy 

E^ / 2ML\^ 



ds^ = (L^-2A/rL-a2)-2/(2-+3)sin-4/(2-+3)^^^2 

+ ■ 2.^, ( 1 - - 2ML - a2)2/(2-+3) sin^/^^^^^) 9[dR' + d9' 



(A + (M2 + a2) sin' 



<l>(i?,^^) = (L^ - 2ML - a2)2/(2-+3) sin4/{2..+3) ^20) 
2aML sin^ ^ , , 



A 



A + a'^ 8111^9 
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Finally by transforming back to the standard Boyer-Lindquist coordinates using eq.(7), we 
arrive at the Bonnor-type dipole solution in Boyer-Lindquist coordinates given by 

/A + sin^eV 



-(Asin2e)-2/(2-+3)rft2 ^ (Asin2e)2/(2-+3) ^ 



(A + (M2 + a2) sin' Of 
2 



+ ( ) (A sin' m{2-+i)/{2-+3), 

^ A + a' sin' 6 ' ^ 



<l> r,e = - Asm'^^ '/('-+3)^ ^= . ^ 2 • 2. ^0 21 
G A + sm 6 

where we restored the "bare" Newton's constant G in the expression for the BD scalar field 

solution. This solution is static and axisymmetric. It is interesting to note that this solution 

becomes the flat Minkowski spacetime asymptotically as r ^ cxo or for M = and a = 

only as u; ^ 00 which amounts to the general relativity limit. And for finite uj, upon setting 

M = = a as well as A = 0, this solution reduces to the empty spacetime solution of 

BD theory given in eq.(9) as it should. Besides, this solution possesses singularities along 



the symmetry axis 9 = 0, vr as well as at r = r+ = M + VM2 + a' where A = for 
finite u. As we shall see later on, these are not curvature singularities but coordinate ones 
on which all the curvature invariants remain finite and indeed the symmetry axis of the 
geometry represented by this solution consists of the semi-infinite lines ^ = 0, vr and the 
segment r = r+. And as we shall see in a moment, at each of the "poles" (r = r+, 6 = 0), 
(r = r+, 6 = 7i), lies two extremal oppositely charged black holes of the BD-Maxwell theory. 
In addition, from the asymptotic (r 00) behavior of the magnetic vector potential given 
above, one can deduce that the solution involves the magnetic dipole moment of 2Ma. Thus 
the solution can be identified with a magnetic dipole solution. Associated with this last 
point, note that changing the sign of a amounts to reversing the orientation of the dipole, so 
we will consider, without any loss of generality, a > 0. Another point worthy of note is that 
even if we take the limit a —>■ 0, this Bonnor-type magnetic dipole solution fails to reduce to 
the single Brans-Dicke-Schwarzschild black hole solution. Noticing that the parameter, a, 
can be regarded (for large a) as indicating the proper separation between the two opposite 
charges |lO, ll^; this observation reveals a peculiar feature that as a — 0, the pair of 
opposite magnetic charges (or black holes) never appear to merge and the limiting solution 
maintains axisymmetric structure. This point is reminiscent of the magnetic black dihole 
solutions in other Einstein-Maxwell-dilaton theories like Kaluza-Klein or low energy string 
theories as well as the general relativity coupled to the Maxwell theory jisj ]. 
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Originally, Bonnor's magnetic dipole solution in Einstein-Maxwell theory was thought to 
describe a singular, point-like dipole. Recently, however, Emparan ^ demonstrated that 
it actually describes a "black dihole" which has regular horizons. Emparan's argument 
for the black dihole interpretation of original Bonnor's solution goes as follows. Bonnor's 
solution is clearly asymptotically flat as r ^ oo and in this asymptotic region, the gauge field 
solution becomes that of a magnetic dipole. Next, although the solution exhibits apparent 
singularities at r = r+ (where A = — 2Mr — = 0), they need to be treated more 
carefully. Note first that the axis of symmetry of the solution (namely the fixed point set of 
the Killing field {d/dcj))) consists of the semi-infinite lines, ^ = 0, vr (running from r = r+ to 
r = oo) and the segment r = r+ (running from 9 = Q io 9 = -k). Then a crucial feature of 
the Bonnor's solution is that at each of the poles, (r = r_|_, 6' = 0) and (r = r+, 9 = n), lies 
a distorted extremal Reissner-Nordstrom black holes. 

Thus it is of equal interest to study if the Bonnor-type dipole solution in BD-Maxwell 
theory constructed in the present work can describe dihole geometry as well and if so, under 
what condition. Among others, the bottomline condition for possible dihole interpretation 
concerns the particular values of the BD u parameter. Note that some time ago it has been 
pointed out [23] that in order for the BD-Reissner-Nordstrom solution (or more generally, 
for the BD-Kerr-Newman solution) to describe a non-trivial, regular, charged black hole 
spacetime, the a;-parameter of the BD theory should take values in the range, —5/2 < uj < 
—3/2. Here, by "regular" it means that the black hole spacetime possesses regular Killing 
horizons at which the invariant curvature polynomials such as R, R^^R^^ , and R^jjafsR^^""^ 
remain nonsingular. Thus hereafter we shall assume that the BD uj parameter takes values 
in this particular range. 

Again, our eventual objective is the dihole interpretation of the BD-Bonnor solution in 
eq.(21). Prior to this, however, it seems necessary to discuss the peculiar geometric structure 
of this BD-Bonnor solution that appears to inherit essentially the same generic features of 
the BD gravity empty spacetime solution we studied earlier in the previous subsection. 
Firstly for finite uj value, this BD-Bonnor solution in eq.(21) fails to be asymptotically flat 
for the following reason. Namely, the solution-generating algorithm of Singh and Rai breaks 
the asymptotic flatness along the way essentially due to the emergence of non-trivial BD 
scalar field representing the inverse of the spacetime-dependent effective Newton's constant. 
This is a generic nature of the Singh and Rai's solution-generating algorithm that we first 
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realized in the construction of the BD gravity empty spacetime solution discussed in the 
previous subsection. Next, the symmetry axis again turns out to be some kind of internal 
infinity. To see this, consider the BD-Bonnor metric solution for the value u — —2 which 
lies in the range, — 5/2 < a; < — 3/2, 

gtt ~ -(Asin^ Qf, Qrr = gee/^ ~ (Asin^ Oy ' ^''"^ ^ (Asin^ 6)^, 
$(r,^) = -^(Asin^^)-^. (22) 

From this, we can readily realize that the symmetry axis is infinite proper distance away 

° y/gfede - /° -^de ^ oo (23) 
Je sm 9 

while the asymptotic region r = oo is finite proper distance away 

/•oo POO 2 

Jr Jr A sm y 

Like before, this can be attributed to the spacetime-dependent behavior of the effective 
Newton's constant 

^-\r, 0) = G'e//(r, e) = ^(Asin^ df. (25) 

Namely, since the Newton's constant is a measure of the strength of the gravitational inter- 
action, this spacetime-dependent behavior of Ge// implies that the gravitational interaction 
vanishes on the symmetry axis while it diverges asymptotically as r ^ oo. More specifi- 
cally, the region such as the symmetry axis is essentially inaccessible as the gravitational 
interaction is effectively absent there whereas the coordinate infinity r = oo can be in a 
finite proper distance as the gravitational interaction grows as ~ and eventually diverges 
at r = oo. That is, the physical (proper) distance between the two spacetime points gets 
larger as the effective gravitational interaction there gets weaker whereas it gets smaller as 
the effective interaction there gets stronger. Again, this is indeed the same generic feature 
of any spacetime solution of the BD gravity theory that wc first encountered in the case of 
the BD gravity empty spacetime solution before. Then next, wc arc led to anticipate that 
the internal infinity nature of the symmetry axis we just observed would likewise lead to the 
occurrence of the conical singularity in the geometry described by this BD-Bonnor solution. 
Thus we turn now to the study of possible conical singularity structure of this BD-Bonnor 
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solution. First observe that the rotational Killing field V'^ = {d/d(j))^ possesses vanishing 
norm, i.e., = gaf3^''^'^ = = [s/(A + sin^ ^)] ^ (A sin^ ^)(2^+i)/(2..+3) = g at the 

locus of r = r+ as well as along the semi-infinite lines ^ = 0, vr for —5/2 < uo < —3/2. This 
implies that r = r+ can be thought of as a part of the symmetry axis of the solution. Namely 
unlike the other familiar axisymmetric solutions, for the case of the BD-Bonnor solution un- 
der consideration (and of course for the original Bonnor solution as well), the endpoints of 
the two semi-axes 6' = and 9 = -k do not come to join at a common point. Instead, the axis 
of symmetry is completed by the segment r = r^. And as 9 varies from to vr, one moves 
along the segment from (r = r+,6' = 0) where a black hole is situated to (r = = vr) 
where the other hole with opposite charge is placed. Obviously, then the natural question 
to be addressed is whether or not the conical singularities arise on different portions of the 
symmetry axis. Thus once again we define the notion of conical angle deficit (or excess). 
Consider that ; if C is the proper length of the circumference around the symmetry axis 
and R is its proper radius, then the occurrence of a conical angle deficit (or excess) 5 would 
manifest itself if {dC / dR)\R-^Q = 277 — 5 jl^. We now proceed to evaluate this conical deficit 
(or excess). The conical angle deficit along the axes 6' = 0, tt and along the segment r = r_|_ 
are given respectively by 

0(o,7t) = 2,71 - I — — \e=o,TT = 2tc, {2b) 



A(l>d^ 



S{r=r+) = 271 - I —— \r=r+ = 271 

where, of course, we used the BD-Bonnor metric solution given above in eq.(21) for u- 
parameter of the BD theory taking values in the range, —5/2 < u < —3/2. Just as we 
expected, we encounter the seemingly pathological "27r-angle deficit" again. This, of course, 
is no surprise as it is not a generic pathology but just another manifestation of the emergence 
of the internal infinity, i.e., the symmetry axis. Namely, as the symmetry axis is infinite 
proper distance away, we, as a consequence, have {gee)^^^^d,^yg^/ d9\g=Q^T, = which, in 
turn, results in 5(o,7r) = Stt. Once again, this conical singularity structure is not a separate 
pathology but a generic feature that appears to inherit essentially the same thing of the BD 
gravity empty spacetime solution we studied earlier in the previous subsection. 
Note that these peculiar features are the generic characteristics of the static solutions of 
the BD-Maxwell theory that appear as long as we confine our interest to the case of the 
black "dihole" interpretation, namely, as long as we assume —5/2 < uj < —3/2 and are 
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really somethings that do not show up in the case of the Einstein-Maxwell theory or other 
scalar-tensor theories such as the minimal Kaluza-Klein theory or the bosonic sector of the 
low energy effective string theories 



Lastly, therefore, we need to ask, prior to its black dihole interpretation, whether or not the 
BD-Bonnor solution given in eq.(21) is qualified at all to be endowed with any meaning as a 
physical spacetime solution. To answer this, recall that for a value of the BD uo parameter in 
the range of interest —5/2 < u; < —3/2, the BD gravity empty spacetime solution was given 
an acceptable physical meaning or interpretation in terms of the spacetime-dependent be- 
havior of the effective Newton's constant. Next, since we keep employing the same solution- 
generating algorithm of Singh and Rai, the Bonnor-type solution of BD-Maxwell theory can 
be thought of as the one that results when one puts a pair of dipole (or oppositely-charged 
black holes as we shall demonstrate shortly) on opposite sides along the symmetry axis of 
the empty spacetime. As can be expected, however, the generic peculiar features of the 
spacetime such as the internal infinity nature of the symmetry axis and its consequence, 
namely the conical singularity still survive the addition process of the dipole content. Thus 
the BD-Bonnor solution deserves essentially the same physical meaning or interpretation as 
that of the BD gravity empty spacetime solution we studied earlier. The only delicate issue 
to be addressed would be whether or not the dipole lying along the symmetry axis is really 
physically meaningful when the symmetry axis is an internal infinity. At least it seems to us 
that it should nevertheless be of enough physical interest as long as we take the BD gravity 
empty spacetime solution seriously as the one that embodies some spirit of Mach's principle, 
namely that the Newton's constant may take different values in different spacetime points. 
This issue also can be viewed from a different angle. Here we are worried about whether or 
not a pair of black holes lying along an infinitely far symmetry axis can be of any meaning. 
Along this line we need to recall an essentially analogous (single) black hole case in which 
the horizon (namely, the black hole itself) is infinite proper distance away to an asymptotic 
observer. In this sense, the BD-Bonnor spacetime is indeed nearly as weird as the familiar 
black hole solutions. 

Since we accepted the BD-Bonnor solution as having a physical meaning, we now proceed to 
elaborate on the magnetic black dihole spacetime interpretation of this solution in some more 
detail. Particulary, we would like to show explicitly that at each of the poles, (r = r_|_, ^ = 0) 
and (r = r+, = tt), lies a distorted extremal BD-Reissner-Nordstrom black holes in a man- 
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ner similar to what happens in the Bonnor's magnetic dihole solution in Einstein-Maxwell 
theory first demonstrated by Emparan. To this end, we begin by performing the change of 
coordinates from {r,6) to {p,6) given by the following transformation law 10, 12, 3| 



r+ + |(1 + cose), smH = ^7=^(1 - cos^^ (27) 



where A(r+) = 0, namely r+ = M + a/ + as introduced earlier. Then one can realize 
that upon changing the coordinates as given above and taking p to be much smaller than 
any other length scale involved so as to get near each pole, the Bonnor-type magnetic dipole 
solution in BD-Maxwell theory given earlier becomes 



-(p2 sin^ ^)-2/(2-+3) M ' ^ (^2 gi^2 ^)2/(2.+3) M ' (^^2 ^ p2^g2 



<l>{p,e) ^ i(p2sin2^12/(2.+3)^ (28) 

Gr 

A ^ q^^==g-\e){l-cose)d<P 

where q = ^^J^^^ and g{6) = [cos^(|) + ^-^^.^"^^2) sin^(|)]. Certainly, this limit of the solution 
can be identified with the "near-horizon limit" (r r+ or p ^ 0) of a distorted extremal 
charged black hole in BD-Maxwell theory with the distortion being described by the factor 
g{6). And due to this distortion factor g{6), the geometry is not spherically symmetric, 
rather it is elongated along the axis in a prolate shape. And of course this deformation of 
the horizon geometry (to a prolate spheroid) is due to the field created by the other hole 
and possibly by the conical defect to which we shall come back in a moment. 
The limiting geometry above were valid for arbitrary value of "a" as long as we remain close 
enough to (one of the) poles. If instead we consider the limit of very large "a" , while keeping 
(r — r_|_) and asin^6' finite, the Bonnor-type magnetic dipole solution reduces, in this time. 



to 



\ -2 



^ - h + ^ j (p2 Sin^ ^")-2/(2-+3)rft2 



+ 1+^ 



2 



(p2 Sin^ 0-)2/(2-+3) ^dp2 ^ p2^Q2^ ^ ^^2 ^^^2 ^)-2/(2c.+3)^2 ^^^2 , 

V PJ ^ 

^pj) ~ i(p2sin2^)2/(2-+3)^ (29) 

Gr 

A ~ g(l — cos9)d(f) 
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with q 



Mr4 



VM2 



M as a — s> oo. 



As we shall demonstrate below, this limit of the solution can be recognized as representing 
the extremal BD-Reissner- Nordstrom black hole (i.e., the charged single black hole solution 
in BD-Maxwell theory). Indeed, this result was rather expected since physically, taking the 
limit a oo, amounts to pushing one of the poles to a large distance and studying the 
geometry of the remaining pole. We now demonstrate in an explicit manner that the above 
limit of the BD-Bonnor solution indeed represents the BD-Reissner-Nordstrom black hole 
spacetime. To this end, we consider the change of the radial coordinate (which takes p to 
the standard Boyer-Lindquist radial coordinate R) given by 



R={p + q). 



(30) 



Then one ends up with 
ds^ -- 



-2/(2w+3) 



RJ 



+ 



R^il 



2/(2a;+3) r 



- 1 



R 



dR^ + R^dtf' 



(31) 



Certainly, this can be identified with the extremal, charged non-rotating black hole solution 
in BD-Maxwell theory found in 



23| . To check this, one only needs to take the non-rotating 
limit of the BD-Kerr-Newman black hole solution there and set r ^ R, 6 ^ 6 and e ^ q 
(see the appendix). Then one gets the BD-Reissner-Nordstrom black hole solution which 
is, in our present coordinates (t, R, 9, 0), given by eq.(31). 

To summarize, based on these observations, we may conclude that the BD-Bonnor solution 
given above indeed describes a black dihole geometry. 



III. Bonnor-type magnetic dipole solution with Melvin-type universe content 

in Brans-Dicke-Maxwell theory 

In the preceding section, we have constructed the Bonnor-type magnetic dipole solution 
in BD-Maxwell theory and demonstrated that indeed it can be identified with a "black di- 
hole" solution particularly for values of the BD cj-parameter in the range —5/2 < uj < —3/2 
representing a configuration in which a pair of oppositely charged extreme regular black 
holes is placed at each of the poles, (r = r^,9 = 0) and (r = r^,9 = n) respectively. 
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However, since we failed to address the detailed nature of conical singularities along the 
symmetry axis, 6 = 0,tt and along the segment, r = for the reason stated in detail 
earlier, we are still away from the satisfactory understanding of how the two oppositely 
charged black holes can be sustained in an (un) stable equilibrium against the combined 
gravitational (same mass) and gauge or electromagnetic (opposite charge) attractions. 
Usually, if there are finite conical angle deficits along the semi-infinite axes, 6 = 0, tt (as 
is the case with generalized Bonnor solutions in Einstein-Maxwell-dilaton theories |l3[), 
one interprets these angle deficits as the presence of cosmic strings providing tension that 
pulls the pole and the antipole at the endpoints apart. Namely, in this interpretation, 
the tension generated by the cosmic strings counterbalances the combined gravitational 
and gauge attractions and holds the dipole (or dihole) apart aginst the collapse. Indeed, 
the emergence of the conical singularities signals the instability of the configuration and 
by the recourse to the presence of cosmic strings with proper tension, one achieves the 
stability as the line r = r+, < < vr joining the poles can be made to be completely 
non-singular then. Alternatively, perhaps another very relevant approach toward the 
stabilization of the system one can turn to would be to introduce an external magnetic 
field, aligned with the axis joining the dipole, to counterbalance the combined gravitational 
and gauge attractive forces by pulling them apart. By properly "tuning" the strength of 
the magnetic field, the internal stresses along the axis would be rendered to vanish. And 
this can be checked in an explicit manner if the introduction of the magnetic field with 
"right" strength leads to the vanishing of the finite conical singularities along the symmetry 
axis, 9 = 0,71 and along the segment, r = r+ Q] at the same time. Particularly 
for the case at hand in which the cosmic string interpretation is not available for the 
reason stated earlier, this second option involving the introduction of external magnetic 
field looks quite appropriate. Note also that the conical singularity structure of the 
generalized Bonnor solution in Einstein-Maxwell-dilaton theories and its cure via the 
introduction of external magnetic field is reminiscent of the Ernst's prescription jJI for the 
elimination of conical singularities of the charged C-metric in Einstein-Maxwell theory. In 
this section, therefore, we attempt to construct the Bonnor-type magnetic dipole solution 
of Brans-Dicke-Maxwell theory containing the magnetic field content that asymptotes to 
the "Melvin-type" magnetic universe which is a flux tube that provides the best possible 
approximation to a self-gravitating uniform magnetic field. And to this end, we shall employ 
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the Ehlers-Harrison-type transformation generalized to the BD-Maxwell theory which, in 
the Einstein- Maxwell theory, has been known to generate an axisymmetric solution with 
a magnetic field content from a known axisymmetric solution without. It, however, is 
interesting to note that there indeed are two different ways of achieving this. The first route 
is to start with the Bonnor-type magnetic dipole solution obtained in the preceding section 
and apply to it the generalized Ehlers-Harrison transformation to arrive at our destination, 
i.e., the BD-Bonnor solution with magnetic Melvin-type universe content. Meanwhile, the 
second route is to start instead with the Bonnor solution with magnetic Melvin universe 
content in Einstein-Maxwell theory (which is already known in the literature [3]) obtained 
from the original Bonnor solution via the Ehlers-Harrisn transformation and then apply 
to it the Singh and Rai's algorithm which is known as a solution-generating method that 
allows one to generate a stationary axisymmetric charged solutions in BD-Maxwell theory 
from their counterparts in Einstein-Maxwell theory. As we shall see in a moment, these 
two methods yield the same final solution and this indicates, among others, that the two 
separate solution-generating algorithms, the Singh and Rai's one and the Ehlers-Harrison 
transformation indeed commute. And in the following we shall discuss both two routes. 



(1) Route 1 



In this subsection, therefore, we construct the BD-Bonnor solution with magnetic Melvin- 
type universe content via the first route stated above. Consider the BD gravity theory 
coupled to Maxwell field described by the action 



S 



1 JcDR ^n^^-n'^Pp F 



We now consider the Weyl-rescaling given by 



(32) 



with n'^{x) 



^0 



{2u + 3) 



(33) 



where Mpi denotes the Planck mass. Under this Weyl-rescaling, the action above transforms 
to 



S = d xJ g 



(34) 
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Namely, upon this Weyl-rescahng, the action of the BD gravity coupled to the Maxwell field 
takes the form of that of the Einstein-Maxwell-dilaton theory in which the dilaton field is a 
minimally-coupled, massless (Weyl-rescaled) BD scalar field. 

Then we are now ready to perform the "generalized" Ehlers-Harrison transformation (already 
known in the literature ^J) which takes an axisymmetric solution of the Einstein-Maxwell- 
dilaton theory to another such solution containing particularly the (asymptotically) Melvin- 
type magnetic universe content. It is given by 

9ij = >^^9i3 ihJT^^), 96^ = ^'^9h^ (35) 



with A = [1 + -BA^]^ + -B^g^<p. 

Indeed, the virtue of the above Weyl-rescaling on the action of BD-Maxwell theory was to 
make use of the generalized Ehlers-Harrison transformation already known in the literature. 
Since we achieved this goal, now the remaining task is to translate this generalized Ehlers- 
Harrison transformation given in terms of the Weyl-rescaled fields {g^u, A^, \E') back into 
that given in the original fields [g^u, A^^, $), which now reads 

9'ij = ^"^gij 9'4>^ = >^'^gH^ (36) 

with A = |1 + ifi^J' + jB'-^S**- 

Finally, upon performing this generalized Ehlers-Harrison transformation (22) on the BD- 
Bonnor solution given earlier in eq.(ll), one gets the BD-Bonnor solution with magnetic 
Melvin-type universe content which is given by 



-(A sin^ ey'^/^'^^+^Ut^ + (A sin^ ^)2/(2c^+3) ^ 



+ A~2(Asin^ 



<(2a;+l)/(2a;+3) 



$(r,^) = ^(Asin2^)2/(2-+3)^ (37) 

Gr 



A = ^— {4Mra + S[(r2-a2)2 + Aa2sin2^]} 



2AS 

.^^ ^ /A + a2sin2^\ 
with A = A 



^ |(A + sin^ 0) + 2BaMr sin^ 9 + ^B^ sin^ ^[(r^ - a'^ f + Aa^ sin^ 9] 
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(2) Route 2 



Next, we attempt to construct the BD-Bonnor solution with magnetic Melvin-type uni- 
verse content via the second route described before. And to this end, we start with the Bon- 
nor solution with magnetic Melvin universe content in Einstein-Maxwell theory discussed 
in and then apply to it the Singh and Rai's algorithm. Again in Boyer-Lndquist-type 
coordinates, the Bonnor solution involving Melvin's magnetic universe content is given by 



'"^^^ + (A + (M2 + a2) sin2 ^^)3 I A + "^^^ 



+ A-^Asin^ 



ds^ = A' 

A = {4Mra + B[{r^ - a^f + Aa^ sin^ 9] } (38) 

with A = ^ |(A + sin^ 9) + 2BaMr sin^ 9 + ^B^ sin^ 9[{r^ - a'^f + Aa^ sin^ ^]| . 



From this point on, the procedure to generate the BD-Bonnor solution with magnetic Melvin- 
type universe content by applying the Singh and Rai's algorithm to this solution is quite 
straightforward and essentially the same as that to generate the BD-Bonnor solution from 
the Bonnor solution in Einstein-Maxwell theory we discussed earlier in section II. To briefly 
sketch the procedure, consider, again, the coordinate transformation given in eq.(17) which 
gives dr'^/A = dR^. Then the Bonnor solution with magnetic Melvin universe content above 
can be cast into the form in eq.(3) from which we can read off 

g2C/s = A^, We = 0, (39) 

yi 

e^^^ = —: ^A^ hl = Asm'^9. 

(A + (M2 + a2) sin^ 9)^ ' ^ 

Then using the solution-generating rule in eq.(5) in the algorithm by Tiwari and Nayak, and 
Singh and Rai, and then transforming back to the standard Boyer-Lindquist coordinates 
using eq.(17), finally we can construct the BD-Bonnor solution with magnetic Melvin-type 
universe content which turns out to be precisely the same as the one obtained earlier via 
the "route 1", namely eq.(37). 

Now let us appreciate the meaning of this agreement once again. If we recall the nature 
of these two methods, the "route 1" involves Singh and Rai's solution-generating method 
followed by the Ehlers- Harrison transformation while the "route 2" involves actions of the 
reversed order. This indicates, among others, that the two separate solution-generating 
methods, that of Singh and Rai and that of Ehlers- Harrison indeed commute. Although 
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this realization may not be so surprising, at least the lesson we learned from it is the fact 
that the Ehlers-Harrison-type transformation we derived in eq.(36) in the context of the 
BD-Maxwell theory was indeed correct and hence can be applied to other cases of interest. 
Now upon successfully constructing the BD-Bonnor solution with magnetic Melvin-type 
universe content, we now turn to the question asking if this solution can describe dihole 
geometry as well. In order to address this issue, we, as before, perform the change of 
coordinates from {r,9) to {p,9) given in eq.(27) and examine the solution for very small 
value of p to get 



„2 „,„2 5N-2/(2a,+3) ( p\ ^^2 ^ (^2 g-^2 ^)2/(2a;+3) ( <^ 



{dp'' + p^dP) 



+ r'(0){p' Sin' e)-'/^^^+'n^] p'sin'Mct)', 



~ i.(^2gi^2^)2/(2c.+3)^ 
Gr 



(40) 



A 



Q 



VM2 + a2 



+ Bq 



g-\e){l- cos e)d(t) 



where q ^ and 0) ^ [cos\^ + (^^ + Bqf sin^ 

charge of the solution is 



+ Bq 



2n 



And lastly, the magnetic 



(41) 



VM2 + a2 

where Acf) denotes the yet unspecified period of the azimuthal angle 0. It is the "physical 
charge" of a single pole computed using Gauss' law Qm = Xsa -F[2] with S' being any 
topological sphere surrounding each pole and F[2] = dA. Again, this limit of the solution 
can be identified with the "near-horizon limit" (r — > r+ or p — > 0) of a distorted extremal 
charged black hole in BD-Maxwell theory with the distortion now being described by the 
factor g{9). And the horizon, i.e., the p — surface lacks spherical symmetry, instead, it is 
a prolate spheroid. 

Lastly, we turn to the conical singularity structure of this BD-Bonnor solution with mag- 
netic Melvin-type universe content. Using the metric solution given above in eq.(37) for 
a;-parameter of the BD theory taking values in the range, —5/2 < cu < —3/2, the conical 
angle deficit along the axes 9 — 0,7r and along the segment r = r+ are given respectively by 

A(j)d, 



27r- 



geedO 



\e=o,n — 27r, 



(42) 
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()(r=r+) = iTT - | , \r=r+ = 27r. 

To our dismay, we are still left with the same embarrassing results as before when the 
magnetic field content was absent. Once again, this seemingly pathological result, namely 
the "27r-angle deficits" along both the semi-infinite lines = 0, vr and the segment r = r+ 
can be attributed to the fact that the metric structure of eq. (37) leading to the emergence of 
the "coordinate singularities" along ^ = 0, tt and r = r+ essentially remained the same (even 
if we have introduced the external magnetic field content into the BD-Bonnor solution) as 
the one of eq.(21) which, as stressed several times, exhibits the internal infinity nature of 
the symmetric axis. Earlier, we mentioned that particularly for the case at hand in which 
the cosmic string interpretation is not available, perhaps this second option involving the 
introduction of external magnetic field could be quite an appropriate approach toward the 
stabilization of the dihole system. And this can be checked in an explicit manner if the 
introduction of the magnetic field with proper strength leads to the vanishing of the conical 
singularities along the symmetry axis, ^ = 0, tt and along the segment, r = r+ at the same 
time. Apparently, however, this hope of ours failed again. As a result, any attempt to remove 
the conical singularities along the symmetry axis and along the segment, r = r+ in terms 
of the cosmic string interpretation or by introducing an external (Melvin-type) magnetic 
field of proper strength failed in the present case of the Bonnor-type spacetime solution 
in BD-Maxwell theory. And we know that this failure can essentially be attributed to the 
internal infinity nature of the symmetry axis to begin with. Nevertheless, we anticipate 
that, although it cannot be checked in an explicit fashion for the reason stated above, the 
introduction of external magnetic field, aligned with the axis joining the dihole, can in 
principle counterbalance the combined gravitational and gauge attractive forces and hence 
eventually stabilize the configuration. 

IV. Summary and Discussions 

In the present work, the construction and extensive analysis of a solution in the context 
of the BD-Maxwell theory representing a pair of static, oppositely-charged extremal black 
holes have been performed. To this end, the algorithm of Singh and Rai's is employed which 
is known to generate stationary, axisymmetric, charged solutions in BD-Maxwell theory 
from the known such solutions in Einstein-Maxwell theory. Indeed it was originally thought 
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that Bonnor's solution in Einstein-Maxwell theory describes a singular point-like magnetic 
dipole. Lately, however, it has been demonstrated [3| that it instead may describe a black 
dihole. Thus it is of equal interest to examine if the Bonnor-type dipole solution in BD- 
Maxwell theory constructed in this work can describe dihole geometry as well. Particularly 
for values of the BD c<j-parameter in the range —5/2 < u < —3/2, it has been demonstrated 
that indeed it can be identified with a black dihole configuration. The supporting argu- 
ment for this conclusion is that although both the BD gravity empty spacetime and the 
BD-Bonnor solution exhibit some peculiar geometric structures such as the internal infinity 
nature of the symmetry axis, they should not be viewed as a pathology but instead should be 
thought of as the manifestation of the Mach's principle that the BD theory itself attempts 
to embody. Then followed the discussion of the nature of this new solution particularly 
concerning its stability issue. Obviously, the configuration described by this black dihole 
solution involves an instability arising from the combined gravitational (same mass) and 
electromagnetic (opposite charge) attractions and it indeed is represented by the emergence 
of conical singularities. As we realized in the text, however, any attempt to address the na- 
ture of possible conical singularity structure of the BD-Bonnor solution along the symmetry 
axis, 6 = 0,TT and along the segment, r = r+ failed essentially due to the internal infinity 
nature of the symmetry axis and the coordinate singularity nature at r = r+ as long as 
we confined our interest to the case of the black dihole interpretation, namely, as long as 
—5/2 < uj < —3/2. As a result, the cosmic string interpretation of the conical singularities 
was not available and therefore we turned to the second option involving the introduction 
of the Melvin's magnetic universe content, i.e., the external magnetic field. Along the way, 
we also noticed that actually there are two possible ways of achieving this and the "route 1" 
involves Singh and Rai's solution-generating method followed by the Ehlers-Harrison trans- 
formation while the "route 2" involves actions of the reversed order. It, then, has been 
realized that the two methods yield the same final solution which indicates that the two 
separate solution-generating methods, that of Singh and Rai and that of Ehlers-Harrison 
indeed commute. Unfortunately, however, this second approach toward the stabilization of 
the dihole system failed again as the metric structure in eq.(37) leading to the emergence of 
the "coordinate singularities" along 6' = 0, vr and r = r+ essentially remained the same even 
if we have introduced the external magnetic field content into the BD-Bonnor solution. This 
failure manifests itself since one cannot determine for the case at hand the proper strength 
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of the magnetic field that can ehminate the conical singularities along the symmetry axis, 
6 = 0,TT and along the segment, r = r+ at the same time. Despite this technical difficulties, 
however, it may be anticipated that in principle the introduction of external magnetic field, 
aligned with the axis joining the dihole, can counterbalance the combined gravitational and 
electromagnetic attractive forces and eventually stabilize the configuration. Lastly, since 
the reinterpretation of Bonnor's solution in Einstein-Maxwell theory as a magnetic black 
dihole configuration, its generalizations in other Einstein-Maxwell-dilaton theories such as 
the Kaluza-Klein or the low-energy string theories have been discussed in the literature 3| ■ 
Thus the only remaining viable scalar-tensor theory of gravity (coupled to Maxwell theory) 
left untouched along this line was the Brans-Dicke-Maxwell theory. And precisely it is this 
issue that we would like to address in the present work. 



Appendix : Regular black hole solutions in Brans-Dicke-Maxwell theory 



In this appendix, we shall briefly review the characteristics of the non-trivial, regular 
black hole solutions of the BD-Maxwell theory constructed and studied in 2^ and particu- 
larly discuss its BD-Reissner- Nordstrom black hole solution limit which is of some relevance 
to the present work. 

The Brans-Dicke-Kerr-Newman (BDKN) solution of the BD-Maxwell theory in Boyer- 
Lindquist coordinates is given by 



(A sin^ 



-2/(2(^+3) 



'A - sin^6'' 



+ 



sm 



2a sin 



A) 



+ (A sin^ 



^ 2/(2(^+3) 



dtdcf) 



G 



(A sin^ 



N2/(2i^+3) 



v4„ 



(43) 



where S = + cos^ 9 and A = — 2Mr + + with M, a, and e denoting the ADM 
mass, angular momentum per unit mass, and the electric charge respectively. Notice that 
the quantities E and A above in eq.(43) are defined differently from those with the same 
names appeared in the BD-Bonnor solution in eq.(21) in the text. Indeed, they are related 
by a — > ia. Just like the way how the BD-Bonnor solution has been derived in the text of 
the present work, this BDKN solution in BD-Maxwell theory has been constructed starting 
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from the Kerr-Newman (KN) solution in Einstein-Maxwell theory and then applying to it 
the Singh and Rai's solution-generating algorithm as well. Note also that the BDKN solution 
above has possible coordinate singularities not only at the outer event horizon where A = 
but also along the symmetry axis = 0, vr. Thus in order to explore the nature of this 
singularity along the symmetry axis, the computation of invariant curvature polynomials is 
necessary. And it is a straightforward matter to realize that the two invariant curvature 
polynomials R^^R^'^ and R^^afiR^^'^^ become finite both on the horizon candidate at which 
A = and along the symmetry axis = 0, vr provided the generic BD tu-parameter takes 
values in the range —5/2 < uj < —3/2 '2^. 

Particularly note that if one takes the non-rotating limit of this BDKN black hole solution 
by setting a (i.e., the angular momentum per unit mass) to zero one ends up with the 
BD-Reissner-Nordstrom black hole solution 
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I sin^^ 
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-1 -2/(2w+3) 



2/(2tj+3) 



- ) dr + sin' 



(44) 



1-'- 

r 



dr^ + r^de^ 



which, upon setting r ^ R, 6 6 and e ^ q, coincides with eq.(31) in the text. 
Now it seems relevant to explore thermodynamics and causal structure of this non-trivial 
BDKN black hole solution in eq.(43) in some more detail. Firstly, these BDKN black 
hole solutions have vanishing surface gravity at the event horizon, k+ = and hence zero 
Hawking temperature, = k^/2tt = provided —5/2 < u < —3/2. In other words, 
they do not radiate and hence are completely "dark and cold". Certainly, this is a very 
bizzare feature in sharp contrast to evaporating black holes in general relativity. Next, 
we turn to their causal structure. The two Killing horizons, i.e., the outer event horizon 
and the inner Cauchy horizon turn out to occur precisely at the same locations (i.e., same 
coordinate distances) as those of KN black hole solutions in Einstein-Maxwell theory, i.e., 
at r± = M ± (M^ 
event horizon at r = r+. 



— e^)^/^. Also it is interesting to note that the proper area of the 



A = [ ded(j){geeg^4>y^^ = ^7r{rl + a^) 



(45) 



is again exactly the same as that of standard KN black hole spacetime. In addition, its 
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angular velocity at the event horizon coincides with that of standard KN solution as well 

- Wshir^) = = -WE\r^). (46) 

Next, observe that the norm of the time translational Killing field 

^"^^ = gu = -{Asm' 9)-'/^'-^''^ 

goes like negative (r_ < r < r+) positive (r+ < r < r^) ^ negative (r > r^) with 
Tg = M + (M^ — cos^ 6* — e^)^/^ > r+ being the larger root of indicating that 

behaves as timelike spacelike timelike correspondingly. And particularly the region 
in which C,'^ stays spacelike extends outside hole's event horizon. This region is the so- 
called "ergoregion" and its outer boundary on which C,'^ becomes null, i.e., r = is called 
"static limit" since inside of which no observer can possibly remain static. Thus if we 
recall the location of the static limit in standard KN black hole solution, we can realize 
that even the locations of ergoregions in two black hole spacetimes, KN and BDKN, are 
the same as well. Namely in two theories, i.e., the BD-Maxwell theory and the Einstein- 
Maxwell theory, rotating, charged black hole solutions turn out to possess identical causal 
structure (i.e., the locations of ring singularities, two Killing horizons and static limits are 
the same) and hence exhibit the same global topology. Thus actually what distinguishes 
the BDKN black hole spacetime from its general relativity's counterpart, i.e., the KN black 
hole is the local geometry alone such as the curvature characterized by the specific a;- values, 
-5/2 <u < -3/2. 

Next, we would like to comment on the divergent behavior of the BD scalar field solution on 
the horizon which, in addition to the null Hawking radiation mentioned earlier, is another 
peculiar feature of the solution. As we mentioned earlier, the BD theory is an alternative 
theory to Einstein gravity and the BD scalar field represents spacetime- varying effective 
Newton's constant, not a matter. Thus the divergent behavior of the BD scalar field in BD 
theory essentially represents the vanishing effective Newton's constant in a certain region 
of spacetime. Besides, since the energy density of the BD scalar field Tj^J^^^^^ vanishes 
and hence satisfies the weak energy condition on the horizon at which A = (of course for 

n 

—5/2 < uj < —3/2) [23], we do not worry too much about the divergent behavior of the BD 
scalar field there. For more detailed account of the nature of the BDKN solution, we refer 



A - sin^ 9 



(47) 



27 



the reader to 



23|. 
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